hep-th/9904197 



Holography, Thermodynamics and Fluctuations of Charged AdS Black Holes 



Andrew Chamblin", Roberto Emparan 6 , Clifford V. Johnson c and Robert C. Myers' 



a D.A.M.T.P., Silver Street, Cambridge, CB3 9EW, UK. 

b Department of Mathematical Sciences, University of Durham, DH1 3LE, UK. 
Departamento de Fisica Teorica, Universidad del Pais Vasco, Apdo. 644, E-48080 Bilbao, Spain. 

c Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506-0055, USA. 

d Physics Department, McGill University, Montreal, PQ, H3A 2T8, Canada. 
Q\ ■ , , 

, a H. A. Chamblin@dampt.cam.ac.uk, Roberto.emparan@durham.ac.uk, c cvj@pa.uky.edu, rcm@hep. physics. mcgill.ca 

Q\ 

5H _ 

Oh, The physical properties of Reissner-Nordstrom black holes in (n+l)-dimensional anti-de Sitter 

spacetime are related, by a holographic map, to the physics of a class of n-dimensional field theories 
coupled to a background global current. Motivated by that fact, and the recent observations of the 
striking similarity between the thermodynamic phase structure of these black holes (in the canonical 
ensemble) and that of the van der Waals-Maxwell liquid-gas system, we explore the physics in more 
detail. We study fluctuations and stability within the equilibrium thermodynamics, examining the 
specific heats and electrical permittivity of the holes, and consider the analogue of the Clayperon 
' equation at the phase boundaries. Consequently, we refine the phase diagrams in the canonical and 

Q\ . grand canonical ensembles. We study the interesting physics in the neighbourhood of the critical 

point in the canonical ensemble. There is a second order phase transition found there, and that 
region is characterized by a Landau-Ginzburg model with A3 potential. The holographically dual 
field theories provide the description of the microscopic degrees of freedom which underlie all of the 
thermodynamics, as can be seen by examining the form of the microscopic fluctuations. 
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I. INTRODUCTION AND SUMMARY 



Explaining the thermodynamic nature of black holes was recognized as an essential hallmark of any complete 
quantum theory of gravity long before such a theory was constructed. The semi-classical approach to quantum 
gravity, which has become quite a mature subject over the years allows for the computation of a number of 

physical quantities. These treatments ignore the details of how a specific solution of Einstein's equations (regarded as 
the effective low energy truncation of the complete quantum gravity) arises, and instead perform a quantum treatment 
of field degrees of freedom in a fixed classical space-time background. 

In that way it was learned that the entropy of Bekenstein H and the temperature of Hawking Q, for example, 
fit into an elegant thermodynamic framework, with questions (such as scattering, unitarity, etc.,) concerning the 
underlying microscopic description — which we might use to construct the underlying "Statistical Mechanics" — best 
left for the future development of a quantum theory of gravity. 

That future is now here. String theory (and/or "M-theory") supplies a microscopic description of the underlying 
degrees of freedom upon which a statistical description of the laws of black hole thermodynamics can be based. This is 
true even though we do not yet have a satisfactory way of writing the theory in all regimes: the "D-brane calculus" (|] 
provides a robust framework within which to describe many properties of black holes || , while in turn being firmly 
rooted in the dynamical framework of string duality and, ultimately, M-theory [Q. Typically, the description of black 
holes (and other important geometrical backgrounds) proceeds by translating their properties into properties of an 
auxiliary field theory, identified as residing on the world-volume of some collection of (D- or M-) branes. 

One of the succinct ways of organizing this microscopic description of the properties of black holes is via "AdS 
Holography" [@H12fl. Then, the thermodynamic properties of black holes in anti-de Sitter space-time are dual to 
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those of a field theory in one dimension fewer jlO],[ll| . The fact that the thermodynamic properties of the AdS black 
holes Jll| are organized by an effective field theory is not implausible, in the light of the fact that AdS acts like a 
natural "box" (with reflecting walls) which neutralizes the tendency of gravitational interactions to render a canonical 
thermodynamic ensemble unstable. The fact that the effective field theory is one which does not contain gravity, and 
that it is actually a "holographically" dual four dimensional gauge theory (with suitable generalizations beyond D—A) 
is another striking example of the fundamental role that gauge theory plays in duality in various situations^. 

That the AdS arena stabilizes the thermodynamics of black holes is especially apparent when one discovers phase 
structures completely analogous to familiar thermodynamic systems from elsewhere in Nature. Such an example can 
be found in the Reissner-Nordstrom-anti-de-Sitter (RNadS) systems in various dimensions ]l7| ]. There, the (Q,T) 
diagram showing the thermally stable phases for a fixed charge (canonical) ensemble turns out to be completely 
analogous to that of the (T, P) phase diagram of the liquid-gas system. The structure of the first order phase 
transitions, etc., are controlled by a "cusp" catastrophe Jig] , common in the theory of discontinuous transitions in 
thermodynamics and many other fields^]. Meanwhile the free energy as a function of temperature F(T), displays the 
characteristic "swallowtail" shape. 

In this paper, we report the results of our further examination of these structures, exploring the equilibrium 
thermodynamics more closely, including the effects of considering electrical stability, and thermal fluctuations. The 
similarities noted between the RNadS physics and that of well-defined systems such as the liquid-gas are more 
than mere analogies: We find that everything has a very natural place in classic equilibrium thermodynamics, as is 
consistent with a holographic duality to thermal field theory without gravity. Accordingly, using the techniques of 
equilibrium thermodynamics, we refine the phase diagrams which we found in ref. [ p"7| somewhat, and identify the 
generic physical properties which give rise to the cusp and swallowtail structures. 

As discussed in our previous paper |HJ , the thermodynamics of the Reissner-Nordstrom black holes in the presence 
of negative cosmological constant in various dimensions are pertinent (because of the holographic map) to the ther- 
modynamics of families of field theories found on the common world-volume of collections of large numbers of branes 
(for example M2- and D3-branes), in the situation where a global background current (or its canonical conjugate 
charge) has been switched on and held fixecl^. 

Geometrically this is performed by simply setting the M2- and D3-branes rotating equally in each of the available 
transverse orthogonal two-planes. The higher dimensional angular momentum becomes the Maxwell U(l) charge 
after the Kaluza-Klein reduction on the (now twisted) sphere, which yields the gauged supergravity. Obtaining a 
pure Maxwell term in this way is not possible starting with the M5-brane, and so the seven dimensional Einstein- 
Maxwell-anti-de-Sitter (EMadSy) theory defines at best a close cousin to the field theory found on the M5-branes' 
world- volumes. The dual theory relevant to EMadSg should be considered in a similar manner. 

A truly rich phase structure for the field theories (with transition temperatures away from T=0) is obtained only 
for finite volume, which is the case we concentrate on here. Our studies correspond to the study of black holes with 
spherical horizons, S n ~ 1 . The field theory resides on IRxS" -1 . The case of infinite volume corresponds to black holes 
with horizons H n_1 , and to field theory on IR". This is of course the case which comes from taking directly the near 
horizon limit of explicit brane solutions. 

As shown in ref. Jl7| , the results for infinite volume may be easily obtained as a scaling limit of the results of 
finite volume, and so we will not discuss them here. Of course, even though we are in finite volume for much of our 
discussion, the thermodynamic limit is still valid here, because the dual field theory is at large N, and a positive power 
of N measures the number of degrees of freedom in the field theory (for example, N 2 in the case of gauge theory, for 
n— 4 here). 

The structure of the paper is as follows: In section || we recall the charged black solutions of the Einstein-Maxwell- 
anti-de -Sit ter system. We also recall the results of performing the Euclidean section and ensuring its regularity. In 
section III , we translate these results into a statement of about the relation between the thermodynamic variables of 



the black hole system in thermodynamic equilibrium i.e., the "equation of state". In section IV, we define the grand 
canonical and canonical thermodynamic ensembles and compute the associated Gibbs and Hclmholtz thermodynamic 
potentials, contrasting the techniques used (and results obtained) to those of our previous work. In particular, we 
note that we can obtain an intrinsic definition of these quantities in Euclidean quantum gravity, by sidestepping some 



1 See also refs. p^-]lq for discussion of how this extends to relating the physics of linear dilaton backgrounds to theories on 
the world-volumes of NS-branes. 
2 Recently, the cusp catastrophe has appeared again the AdS/CFT context, in ref. jHj]. 
3 See |J(J for additional work on how to relate charged AdS black holes to string/M-theory. 
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of the technical subtleties — encountered in the "background subtraction" technique for regularizing the action — in 
favour of the "counterterm subtraction" technique [|2l|]22[ | . In the rest of the section, we examine the features of these 
potentials quite closely, in preparation for later detailed studies. In section [v|, we use the equation of state and the 
first law of thermodynamics to identify the origins of the crucial features of the shape of the Hclmholtz potential 
(free energy) . This "swallowtail" shape is responsible for the interesting phase structure in the canonical ensemble. 



Section |V| examines the conditions for thermodynamic stability of the black holes, examining the specific heats and 
permittivity of the black holes. In this way, we identify the stable regions of the solution space of the equation of state. 
We use this stability information, together with the information gained in earlier sections, to deduce the refined phase 



diagrams exhibited in section VII , and some details of the phas e dia grams (the slope and convexity of the coexistence 



curves) are refined by using the Clayperon equation in section VIII . 

As already stressed in this section, the thermodynamic quantities and studies performed in those sections are rooted 
firmly in a microscopic description. This is ensured by the fact that we can in principle embed this entire discussion 
into a complete theory of quantum gravity: string (and/or) M-theory. In practical terms, this microscopic description 
— the "statistical mechanics" underlying the thermodynamics — is summarized neatly in terms of the holographically 
dual field theory. In this way, therefore we may carry our calculations further and examine the nature and magnitude 
of the microscopic fluctuations of the various thermodynamic quantities we have computed, knowing that we have a 



description of their origin in field theory. Thus, we find in section IX that the fluctuations behave in a way consistent 
with the underlying microscopic physics being supplied by the field theory: the size of the (squared) fluctuations is 
controlled by a prefactor which corresponds to precisely the inverse of the number of degrees of freedom of the dual 
field theory. We observe that the size of the fluctuations diverges as the system approaches a critical point in the 
(Q, T) plane. 

Through most of the paper, we carried out our computations for the four dimensional case, in order to keep many 
of our formulae simple. Section [x| collects together some of the results for the computation of various quantities. We 
stress that the qualitative structure of the physics is the same for all dimensions d>4, where d—n+1. Briefly, we 
also discuss in that same section the issue of the meaning of the formal definition of other thermodynamic ensembles 
by Legendre transform. It is not always the case that the thermodynamic quantities thus defined may be arrived at 
by (known) computations in Euclidean quantum gravity. Therefore, interpretations of the physics of such ensembles 
are to be taken with (at least) a pinch of salt, until such time as new technology becomes available to compute the 
relevant quantities directly in quantum gravity, as we have done here for the fixed potential (grand canonical) and 
fixed charge (canonical) ensembles. 



Section XI discusses the underlying structure of the phase structure of the canonical ensemble in the neighbourhood 
of the critical point. In particular, the physics local to critical point is universal for all of dimensions d>4. The critical 
point is a second order phase transition point at the end of a coexistence line of first order phase transitions. As such, 
it has a universal description in terms of a Landau-Ginzburg model, with a quartic potential — A3 in the A-D—E 
classification of such potentials. The deformation of this potential gives the classic "cusp" catastrophe which underlies 
the critical behaviour, as is well known from the van-der Waals-Maxwell descri ptio n of the liquid-gas system, with 
which our black hole physics shares many features, as originally reported in ref. |17j]. 

In closing the introduction, we would like to stress once again how elegantly the properties of anti-de Sitter space 
yield charged black hole physics so closely akin in structure to that of ordinary field theory-like systems, with which 
we have more intuition. 

From the point of view of the Maxwell part of the action, the black holes are nothing more than spherical capacitors, 
and as such, the amount of energy they can store grows with the charge on them, but falls with increasing hole radius. 
From the point of view of the Einstein-Hubert action however, the black holes store an amount of energy which grows 
with radius. After a little thought, one might expect on general grounds, therefore, that there might be an interesting 
phase structure resulting from a competition between these two pieces of the action. 

Such reasoning on its own would not be enough to genuinely fill out the whole {Q, T) phase diagram, as the equation 
of state needs additional structure. It is the presence of (negative) cosmological constant which provides this final 
part: First, it provides black hole solutions which are thermally stable in ensembles involving fixed temperature [fl3|| , 
but secondly, as it defines a new length scale, is allows the system to distinguish, on the one hand, black holes which 
are large from those which are small, and on the other hand, black holes which have small charge, from those with 
large charge. 

It is because of these features that the charged black hole thermodynamics has a chance to be similar to the van-der 
Waals model of the liquid-gas. Recall that without the inclusion of the effects of the length scales set by finite particle 
size, on the one hand, and attractive inter-particle forces on the other, that system would have only the much less 
interesting physics of the ideal gas: there would be no competing effects, as a function of length scale, with which to 
trigger a phase transition. These basic features of AdS give holography a chance to work in a way which is consistent 
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with our intuition that the microscopic physics should be modelled by ordinary field theory. 



II. CHARGED ADS BLACK HOLES 



For space-time dimension n+1, the Einstein-Maxwell-anti-deSitter (EMadS„+i) action may be written asj 

n[n — 1) 



16ttG 



M 



R — F 



r- 



(i) 



with A=— n ^2p being the cosmological constant associated with the characteristic length scale I. Then the metric 
on the Reissner-Nordstrom-anti-deSitter (RNadS) solution may be written in static coordinates as |23| , p^ , p7| 

d«a = -V(r)dt a + ^ + r a d^_ 1 , (2) 

where dCl^—i is the metric on the round unit (n— l)-sphere, and the function V(r) takes the form 

nr) = l-^ + -i^ + ^. (3) 
Here, m is related to the ADM mass of the hole, M (appropriately generalized to geometries asymptotic to AdS p5fl), 



as 



(n- 1)^-1 

where u> n -i is the volume of the unit (n— l)-sphere. The parameter q yields the charge 

Q = y/2(n-l)(n-2)(^)q, (5) 

of the (pure electric) gauge potential, which is 

A =(-- c -h+*)dt, (6) 

where 



P) 

n — 1 

and $ is a constant (to be fixed below). If r + is the largest real positive root of V(r), then in order for this RNadS 
metric to describe a charged black hole with a non-singular horizon at r=r +1 the latter must satisfy 

r 2n-2 + p^n-i > q 2j2 (8) 



71-2 



Finally, we choose 



$ = (9) 



4 We scale the gauge field so as to absorb the prefactors involving the (7(1) gauge coupling into the action. 
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which then fixes A t (r + )—0, as is required by (Euclidean) regularity of the one-form potential (||) at the fixed point 
set of the Killing vector 9 t . The physical significance of the quantity $, which plays an important role later, is that 
it is the electrostatic potential difference between the horizon and infinity. 

If the inequality in eqn. (||) is saturated, the horizon is degenerate and we get an extremal black hole. This inequality 
imposes a bound on the black hole mass parameter of the form m>m e (q, I). 

In passing to the thermodynamic discussion, we define the Euclidean section (t—>ir) of the solution, and identify the 
period, (3, of the imaginary time with the inverse temperature. Using the usual formula for the period, /?=47r/V (r+), 
which arises from the requirement of regularity of the solution, we obtain: 

47rl 2 r 2 + n - a 



nr. 



+ (n - 2)Prf- A - (n - 2)q 2 l 2 



This may be rewritten in terms of the potential as: 

4nl 2 r + 

= (n-2)l 2 {l-c 2 <S> 2 )+nr 2 + ' 

For simplicity, we will specialize to n=3 (therefore working with EMadS.4) to avoid cluttering our expressions with 
complicated dependences on n. Our results will remain qualitatively the same for higher n (see the comments in 
section XI ) , and we list some of the n-dependent formulae in section M. Our analysis is further simplified by adopting 



the following rescalings (once we have set n=3): 



V3 " I 

and for the various thermodynamic quantities used in ref. fll 



T^^T, Q_>^GQ, (12) 



{W,F,E}->^G{W,F,E} , S^-^S. (13) 



and r + — » — — r+ . (14) 

Essentially we are introducing a system of dimensionless quantities in which everything is measured in units of the 
AdS scale I. This scaling is chosen so that the thermodynamic formulae still all have their standard form, i.e., 

dE = TdS + <PdQ . dF = -SdT + <S>dQ dW = -SdT ~ Qd$ , etc. (15) 
In the following, all of the quantities which follow are the rescaled dimensionless quantities, unless stated otherwise. 



III. EQUATION OF STATE 



The Euclidean regularity at the horizon discussed at eqn. (nffl) is equivalent to the condition that the black hole 
is in thermodynamical equilibrium. The resulting equation ( |l0| ) may therefore be written as an equation of state 
T=T($, Q) (analogous to T=T(P, V) for, say, a gas at pressure P and volume V). For n—3, one finds: 

$ 2 (1 _ $2) + Q 2 

2Q$ ' ^ ' 

One can also solve for Q as 



Q = T$ ± $ VT 2 + $ 2 - 1 . (17) 

From this equation of state we see that for fixed $ we get two branches, one for each sign, when the discriminant under 
the square root is positive. For fixed Q, T($) has three branches for Q<Q C rit and one for Q>Q C rit, where the critical 
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charge is determined solving for the "point of inflection" where (dQ/d$) T — (d 2 Q/d$ 2 ) T =0. In the dimensionless 

units used here, one finds Q cr i t =l/(2v / 3), X , C rit=2v / 2/3, $ crit =l/v6, E CT it=V2/3, and r + ( crit )=l/v / 2- It is useful to 
plot the ($, Q) isotherms, i.e., plot Q( < E>) for fixed T, and we exhibit these in figure |[ 




FIG. 1. Plots of the equation of state of $ vs. Q, showing isotherms above and below the critical temperature T cr it. 
For T<T cr it, there is only one branch of solutions, while for T>T cr i t , there are three branches. The values of T for the 
isotherms plotted are (top down) T=0, 0.8, Tci-it, 1.0, 1.2. The central (dotted) curve is at the critical temperature. 



As T goes to zero, we approach the extremal black holes. Their equation of state is 

§1 = -(1 + v 7 ! + 4Q 2 ), for arbitrary T. (18) 

For some later computations, it is often convenient to use as an additional, non-independent parameter, the black 
hole radius r + , in terms of which 



1 



Q = r + ^/r 2 + -2r+T+l, 2T = r+ + — - ^- , (19) 



and 



* = ^- = v'Vi - 2r+T + 1 . (20) 



IV. THE GRAND CANONICAL AND CANONICAL ENSEMBLES 



In thermodynamic parlance, the "grand canonical ensemble" is defined by coupling the system to energy and charge 
reservoirs at fixed temperature T and potential $ (an intensive variable). The associated thermodynamic potential is 
the Gibbs free energy, W[T, $>]=E—TS—$Q. Holding the extensive variable, Q, fixed, on the other hand, defines the 
canonical ensemble, with its associated thermodynamic potential is the Helmholtz free energy F[T,Q]=E—TS. See 
section [x| for a brief discussion of other ensembles. 

In ref. jL7|, the calculation at fixed potential was carried out by computing the action d la Gibbons-Hawking. With 
that technique, one must regularize the computation (as the action is formally infinite) by subtracting a contribution 
from a "reference" background which matches the solution of interest asymptotically, giving a definition of the action 
relative to that of the reference spacetime. In this case it is appropriate to use AdS — with a fixed (pure gauge) 
potential at infinity — as the reference background. 



G 



Remarkably AdS spacetime provides for another regularization which yields an intrinsic definition of the action. In 
other words, the computation makes no reference to any other solution of the equations of motion. Instead, the method 
pl| , p2t proceeds by adding a series of boundary counterterms to the action. We refer to this as the "counterterm 
subtraction" method of defining the action, a technique tailored to spacetimes which are locally asymptotic to anti-de 
Sitter, as the counterterms are defined on the natural boundary, with which such spaces are endowed, using the AdS 
scale I. Also note that the inclusion of additional sectors to the gravitational and cosmological parts of the action, 
such as Maxwell terms, does not affect the definitions and therefore we can still use the same counterterms in the 
present context. 

The results, using either the reference background or the counterterm subtraction methods, are identical for the 
particular case in which we want to fix the potential^, since it is possible to have AdS space as a background solution 
at arbitrary temperature and (constant) potential (but, crucially, see later). In the present notation, the answer is: 



W[$,T\ = 



1 

12 



3§(1 



$2)_ 



(21) 



Here, Q is given as Q($,T) by the equation of state (|l7|). In terms of r + , this is 

W=^[ir + {l-<5> 2 )-r\] , 



(22) 



and it is plotted in figure ^, with choice slices displayed in figure ^. 

1 




FIG. 2. Plots of the Gibbs potential W[$,T] in three dimensions. 



5 For even values of n there appears a Casimir energy term |21 
here. 



which is immaterial for the discussion of thermodynamics 
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Turning to the canonical (fixed charge) ensemble, we wish to compute the Helmholtz potential F[Q, T] (a.k.a, the 
"free energy"). In ref. Jl7t , where we used the reference background method to compute this, it was necessary to 
compute the action using an extremal black hole as the reference background. This is because anti-de Sitter space 
with a fixed charge Q, as measured at infinity is not a solution of the equations of motion and so is not an appropriate 
background. In order to get an intrinsic definition of the action for fixed charge therefore, we employ the method of 
counterterm subtraction, yielding: 



3 i- If) +w 



(23) 



where $ is given as $(Q,T) by the equation of state (|l7|). In terms of r+, F[Q, T] may be written 



F[Q,T] 



12 



Sr. 



,91 

r+ 



(24) 



As a consistency check that we have performed the computation correctly, note that this result may be obtained 
from the result for the Gibbs potential by formally calculating the Legendre transform F[Q, T]=W[$>, T]+Q&. When 
computing F from a Euclidean action, the additional Q$ term has its origin in the boundary term introduced so 
as to recover the correct variational problem from the action. It is especially satisfying to see that the counterterm 
subtraction method places such intuitive relationships from equilibrium thermodynamics on a firm footing. We shall 
have more to say about this in section |x| 

In ref. |l7j], where we computed the action using an extremal reference background, we obtained the following 
expression for the free energy (which we denote here as F): 



F[Q,T] 



12 



9Q$ - 4^- - 



(25) 



Note that in this case, one should consider <!>=$— <I> e as the state variable, instead of $. Then, the first law is in this 
case dE—TdS+^dQ, and E measures the energy above the extremal state. Furthermore, it is with $ that W[T, $ 1 of 
eqn. (21) and F[T, Q] are Legendre transforms of each other, as they should be. While F[Q, T] as computed in ref. jlj 



using the extremal background is in no way problematic, we shall not examine it further here, as the new technology 
of the counterterm subtraction method has supplied us with an intrinsic definition of the Helmholtz potential, which 
is the more natural Legendre-transform partner of the Gibbs potential ( ^l| ) found earlier. 

We shall see that the qualitative features of the results obtained in ref. jl7| for the canonical ensemble using 
F[Q, T] will persist here, as the extremal background subtraction essentially redefines the absolute normalization of 
some results. (The later analysis of intrinsic stability which we do in section VI would have to be somewhat modified 
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before direct comparison to the extremal subtraction results, however, as we will make heavy use of the equation of 
state in terms of variables (<i>,Q,T), and not the triple (<S>,Q,T) appropriate to that case.) 

We now return to the analysis of the intrinsically defined Helmholtz potential F[Q,T}. It was noticed in ref. J!?]] 
that a plot of F(T) for various values of Q reveals (below a Q C rit) a section of a "swallowtail" shape, which controls 
much of the phase structure (in the canonical ensemble) discussed there, and to be discussed later here. (See figs. 5 
and 6 of ref. jl7| and associated text for details.) The same may be observed here for F(T) for varying Q, as shown 
in figure ^. 



0.4 -i 0.4 -i 0.4 -i 




-0.3 J ' -0.3 J ' -0.3 J 



FIG. 4. The free energy vs. temperature for the fixed charge ensemble, in a series of snapshots for varying charge, for values 
Q=0, 0.15 and Q=0.299. Note that Q cr it=0.289, so in the last plot, the bend (near T cr i t =0.943, is in the neighbourhood of the 
critical point of second order. 

It may be further observed that a plot of F(Q) for fixed T reveals (above a T C rit), a similar swallowtail section, as 
shown in figure |[ 







0.3 0.4 

Q 




-0.2 J _ -0.2 J -0.2 J 

FIG. 5. The free energy vs. charge for the fixed charge ensemble, in a series of snapshots for varying temperature, for values 
T=0. 943, 0.997, 1.00, 1.10, and (for the "Zorro" plot) T=T H p=1.154, and finally T=1.20. Note that T cri t =0.943, and so in the 
first plot, the bend (near Q cr i t =0.289), is in the neighbourhood of the critical point of second order. 



The full three dimensional shape of F[Q,T] is plotted in figure^. 
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FIG. 6. Plots of the Helmholtz potential F[Q,T], in three dimensions, clearly showing the swallowtail shape for T>T CIlt and 

Q<Qcrit- 

That such a shape appears in the thermodynamics (above T cr ;t or below a Q C rit) can be shown to follow from the 
first law of thermodynamics, the definition of the thermodynamic potentials, and the form of the particular equations 
of state which the black holes obey. We will show how this comes about next. 

V. THE SWALLOW TALES 

The sections of swallowtails in the F(Q) and F(T) plots (above a T cr it or below a Qcrit) can be seen to come from 
the existence of the previously mentioned three branches of solutions to the equation of state. We have from the first 
law, and the definition of the thermodynamic potential, that dF=—SdT+$dQ. Therefore, for fixed T we find 



F(T)= J <S>(Q)dQ + f(T) , (26) 

where f(T) is an arbitrary function of T. The integral function can be obtained by looking at the plot of isotherms. 
When we have three branches (i.e., T>T cr i t ), the curve &(Q) winds back and forth in a way that the integral 
describes a shape with three connected branches, constituting a section of the "swallowtail" shape. This can be seen 
by examination of the plots of the equation of state in figure [j] and the plots of the slices F(T) displayed in figure |[| 

Equation ( p6| ) is usually employed to formulate an "equal area law" governing the phase transitions of the system. 
The latter occur at the point where the free energies of two branches, (say A and B), are equal: Fa—Fb- From 



eqn. (26) this equality may be translated into a statement about the equality of the areas enclosed by the isotherm 
curves and a line of constant Q in the ($, Q) plane, as shown in the sketch on the left in figure ^. There is a subtlety, 
though, in using eqn. (^6| ) with the isotherm curves of eqn. ( Ji7| ) for T>1 (recall that isotherms with T>1 go through 
the origin <j>=Q=0. See figure n .). Given that the transition is governed by the equal area law, it would seem from 
the curves on the right in figure 7] and the area law deduced from eqn. (^6|) that even for T>1, for which a minimum 
value of $ ceases to exist, one can always find a phase transition point for arbitrarily large temperatures and small 
enough charge. This must be wrong since it contradicts what we know about the phase transition from the curves of 



6 As visual differentiation is often easier to perform than integration, we gently remind the reader that the defining relation 
&=(dF ' j '8Q)t may be of use here, in conjunction with the snapshots of F(Q) for fixed T given in figure |E| 
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F for constant Q, namely, that the phase transition takes place at a temperature t hat is smaller (or equal, at Q=0) 
than the Hawking-Page temperature Tup. (See ref. Jl7j and the upcoming section VII for detailed discussion of the 
phase structure.) 




Q Q 

FIG. 7. The figure on the left shows how the condition for a phase transition may be interpreted in terms of an "equal area 
law" analogous to that due to Maxwell for the van-der Waals liquid-gas model. For T>1, though, the isotherms have a very 
different qualitative structure. The equal area law one might formulate, deducing a phase transition for arbitrarily high T>1, 
for small enough Q, actually is incorrect. See text for the resolution of the puzzle. 



The resolution of this puzzle is instructive, and is made manifest most clearly by working in terms of the parame- 
ter r + , using eqns. (|l9|), ( p(i|) and (|2~i|). We can explicitly compute F using eqn. (26) as 



F= / ^(f+)(dQ/df + )df + (27) 
Jo 

to recover precisely eqn. (|^). So, what is the reason that "naive integration" using the "the equal area law" yields a 
different result? 

The point is that for T>1 the function F(Q) is discontinuous at Q=0, where branches 2 and 3 separate (see for 
example, the last plot of figure ||). For those isotherms, there is a range of values for r+, T-VT T ^l<r+<T + VT^I 
for which Q and $ become imaginary. Nonetheless, the product QdQ is real throughout, and so is F. Then, F(Q) 
would be a continuous function if we plotted it in the complex Q plane. In performing the integration above for T>1 
we have implicitly included the points where $ and Q are imaginary. Notice that it is by including these points that 
we recover sensible physics, since we want the critical line to end at Q=0 at the point of the Hawking-Page phase 
transition. The "equal area law" , as it is, fails in this instance. 

Let us now turn to the study of the free energy for fixed Q. We have 

F(Q) = - [ S(T)dT + g(Q) (28) 



In this case we need S(T). Since 



1 /<3\ 2 A 



S =2 f "f < 29 > 



we can use the equation of state to plot S(T) for fixed Q, which is shown in figure |[ 
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s 

FIG. 8. Plots of the equation of state of S vs. T, showing isocharge lines above and below the critical charge Q cr it. For 
<2>Qcrit, there is only one branch of solutions, while for Q<Q C rit, there are three branches. The values of Q for the isocharge 
curves plotted are (top down) Q=0,0.20, Q C rit,0.45, 0.80. (The uppermost curve shows the Q=0 case, which has two branches. 
The central (dotted) curve is at the critical charge. 

It can be readily seen that for Q<Q C rit we get three branches (notice that the qualitative features of the plot 
of S(T) follow from those of r+(T) or r+((3) plotted in figure 3 of rcf. ]l7| , where it has a resemblance to the van der 
Waals P{V) curve). A section of the swallowtail again follows^ 

The astute reader may wonder why the swallowtail shape (and the resulting liquid-gas-like) phase diagram occurs 
in the canonical ensemble, where in addition to T, the extensive variable Q is an external control parameter, and not 
in the grand canonical ensemble, where the intensive variable $ would be the control. This is of course what happens 
in the case of the van-der-Waals-Maxwell system, where the phase diagram is in (P, T) space, and not (V, T) space 
p6|p7[ . The swallowtail shapes occur there in the Gibbs potential. It is now hopefully clear that the answer follows 
from the fact that our equation of state yields three branches of solutions for the intensive variable <E> (or T) as a 
function of the fixed extensive variable Q (or S), as can be seen by examining the curves displayed in figs. [I] and |^. 

That there are no swallowtail shapes in any of the other ensembles follows from the fact that no more than two 
branches occur for the equation of state written in terms of other variables. 



VI. INTRINSIC STABILITY 



Given that we have the full power of the thermodynamic framework at our disposal, (thanks to the stabilizing 
influence of a negative cosmological constant), it is interesting to consider the thermodynamic stability of our various 
solutions against microscopic fluctuations^. Notice that one can always formally compute the relevant macroscopic 
quantities (like specific heats, etc.,) which we discuss here, without any reference to an underlying microscopic descrip- 
tion. This has been done in the context of black hole thermodynamics since time immemorial. The difference here 
is that we know the nature of the microscopic degrees of freedom which supply the underlying "Statistical Mechanics " 
which gives rise to these macroscopic thermodynamics quantities. The underlying physics is that of the gauge theory 



7 Again, one can use figure |^ to reconstruct F(Q) visually using the integral relation, or one may use the definition of the 
entropy S= (dF/d T)Q to reconstruct figure |§| from figure ^J. 
8 See refs. pSlp^l for analyses which overlap with those presented here, in a similar context. 
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to which this system is holographically dual, which in turn is the physics of coincident branes. This will become more 



apparent in section [X when we explicitly study the fluctuations themselves. 



Thermodynamic stability may be phrased in many different ways [[30 31|, depending on which thermodynamic 
function we choose to use, and how obscure we are attempting to seem. For example, it can be seen as minimization 
of the energy, E, as a function of (S, Q), or maximization of the entropy S, as a function of (E, Q), etc. In any case, 
one is considering an infinitesimal variation of the state function away from equilibrium. The first law (|l^) will ensure 
that the first order terms vanish. Stability is then a statement about the second order variations. Generally then 
the stability conditions are phrased in terms of the restriction that the Hessian of the state function is positive (or 
negative, depending on the context) semi-definite. 

An equivalent but physically more transparent way of writing the stability conditions is in terms of specific heats 
and other "compressibilities", to wit: 

^m v >-°- HIX^ 

The first two, the specific heats at constant electric charge and potential, are familiar analogues of the specific heats 
at constant volume and pressure in fluid systems. In the case in hand, they determine the thermal stability of the 
black holes, indicating whether a thermal fluctuation results in an increase or decrease in the size of the black hole. 
(This follows from the fact that the entropy is proportional to the size of the black hole,.) Stability follows from C>0, 
given the fact that black holes radiate at higher temperatures when they are smaller. 

The last quantity, st, has the following physical interpretation. It is negative if the black hole is electrically 
unstable to electrical fluctuations (if they are possible, see later discussion). This happens if the potential of the black 
hole decreases as a result of placing more charge on it. The potential should of course increase, in an attempt to 
make it harder to move the system from equilibrium^]. Et therefore deserves to be called the "isothermal (relative) 
permittivity" of the black hole. 

There are of course other interesting "response functions" for the system, such as the adiabatic permittivity, 
(dQ/d$)s or the quantity analogous to the coefficient of thermal expansion in liquid-gas systems, a$=(<9Q/9T)$, 
which are not all independent. The ones which we have discussed above will suffice for the physics that we study in 
this paper. 

We may examine the plots of the isotherms in figure [I] and deduce that the negatively sloped branches are electrically 
unstable if there are electrical fluctuations possible. Similarly, we may deduce that the negatively sloped branches of 
the (S, T) isocharge curves in figure ||, are thermally unstable, and so on. 

Stability follows, equivalently, from the concavity /convexity of the plots of F and W as functions of T. In fact, the 
specific heat conditions are equivalent to 

<« ^d (^) <0, (31) 



d 2 F\ ( 



whereas the permittivity condition is 



d 2 F\ fd 2 W\ 

W )^0, 0, (^) t <0. (32) 



By examining the isotherms displayed in figure ^, we see that there are a number of features in the (Q, T) plane 
which govern electrical stability. Gcncrically, let us describe the three branches of an isotherm as follows: We call 
"branch 3" the branch of solutions which extends all the way from Q=oo, terminating where dQ/d$>=Q. From there, 
"branch 2" takes over, terminating where again dQ/d$>=0. The isotherm continues with "branch 1" until the point 
Q=0, $=1 is reached. This terminology matches that of ref. 

From this definition, then, branch 3 is electrically stable for most of its extent, except for a small region near the 
join with branch 2. In this case, before reaching the point where dQ/d&=0 the permittivity changes sign at a point 
where dQ / d^—oo and renders branch 3 electrically unstable thereafter. This is a feature that is absent from the 



9 This follows from common sense, or more formally, Le Chatelier's principle. 
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standard van-der-Waals-Maxwell system (in the latter there are no points in the isotherms where dP/dV — oo), and 
which will introduce a significant modification of the phase diagram. 

Branch 2, being between two places where dQ/d&=0, has positive definite slope and hence is electrically stable 
everywhere, while branch 1 is electrically unstable everywhere, having negative definite slope. To compute precisely 
where the electrical instability begins, we need only find the location of the minimum of the isotherms, that is, the 
above mentioned point where dQ/d^—oo, which is given by the equation Q=T\J\—T 2 . With the segment of the 
T-axis from to 1, this forms a region in the (Q, T) plane within which branch 3 and branch 1 are unstable to electric 
fluctuations. Branch 2 is electrically stable everywhere, as mentioned before, but as already pointed out in ref. [ jl7| , 
and as a quick examination of the figure || of the isocharge (S, T) curves reveals, branch 2 is unstable to thermal 
fluctuations, and so never plays a role in the canonical and grand canonical ensembles. 

It is also entertaining to subject by eye the snapshots of F and W taken in figures 0, |] and || to the convexity and 
concavity conditions ( J3l| ) and ( p2] ) . We find that the shapes of F and W do indeed confirm our conclusions about the 
stability of the various branches. 

It is very instructive to plot the boundaries of the various branches in the (Q,T) plane: 




1 j 7 T 

FIG. 9. The demarkation of the various branches of black holes in the (Q,T) plane. Points on branches 1 and 3 which 
lie inside the solid curved line are unstable to electric fluctuations. Branch 2 is electrically stable but thermally unstable 
everywhere. 

It is particularly interesting to note that the figures in the previous plot are simply the three sheets of an underlying 
"cusp catastrophe" shape, as can be seen by assembling them in three dimensions to reconstruct the equation of state 
in figure |l|. Indeed, it is highly instructive to align the surface $(Q,T) describing equation of state and the surface 
F[Q, T] giving the swallowtail shape of the free energy, in such a way as to project some of their important features 
to the (Q, T) plane, as done in figure [l^. This gives rise to the critical phase diagram which we will discuss in the 
next section. 
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T 

FIG. 10. The swallowtail shape (free energy) and cusp shape (equation of state) for the charged black hole thermodynamic 
system. Notice the features which result in the critical line and point in the (Q, T) plane. 



As anticipated, the shape formed by the equation of state in the neighbourhood of the critical point is merely a 
distortion of the standard cusp shape, which was encountered in the variables (r+,Q,/?) in our previous paper |L7|. 



Figure 
section 



111 shows this standard shape with two sample trajectories in state space. 



It will be discussed further in 



1G 



FIG. 11. A sketch of the cusp catastrophe in action (in (r+,Q,/3) space). Two sample trajectories are shown, one (Q<Q C rit) 
encountering a phase transition, the other (Q>Q C rit) does not. The precise location of the line across which the transition 
happens is given by the minimum free energy condition, or equivalently, an appropriately formulated "equal area law" . 

As a final comment that in cases where one of the local stability criteria (|30| ) are violated, we are not always able 
to determine the stable ground state. However, the precise nature of the stability violation is providing information 
about how the system will relax to a new stable configuration. For example, one has 5 2 E(xst SQ 2 and so et<0 
indicates that the black hole should relax by reducing its charge, i.e., it will emit charged particles (if possible). 



VII. PHASE STRUCTURE 



Figures [T(] and ^, together with the slices displayed in figures [| and ||, show how the free energy curve determines 
the phase structure of the black holes as one moves around on the state curve in the canonical ensemble, while figures ^ 
and the slices displayed in figure ^| determine the phase diagram for the grand canonical ensemble. We performed 
this analysis in ref. [ Il7j , and we recall it here for completeness, before going on to refine the resulting phase diagram 
using the information uncovered in this paper. 

The dashed line in the (Q,T) plane shows the boundary of the region multiply covered by in the state curve, 
and correspondingly, the free energy has three possible values in that region also (see fig. |10|), which constitutes the 
swallowtail region. The free energy of branch 2 is always greater than that of either branches 1 or 3, however, and 
so there is no transition along the dashed lines. Along the solid line, the free energies of branches 1 and 3 are equal, 
and there is a first order phase transition (the first derivative of the free energy is discontinuous) along this line. 
Also note that the one dimensional Q=0 situation is the familiar Hawking-Page transition ]l3| between AdS and 
AdS-Schwarzschild, which happens (in our units) at T=T#p=2/\/3~1.154, for n=3. 

The solid line is the "coexistence curve" of the two phase of allowed black hole. The line ends in a critical point. 
Above this point, there is no transition, and one goes from large to small black holes continuously (the distinction 
between branch 1 and branch 3 is removed) . (The reader should compare this to the physics of the liquid-gas system 
for an exact analogue in classic thermodynamics.) As the first derivative (but not the second) of the free energy F is 
continuous at the critical point, there is a second order phase transition there, about which we will have some more 
to say in sections 



IX and XI. 



This physics is all summarized in figure |12|, where we have also displayed the phase diagram in the grand canonical 
ensemble (the (3>,T) plane), which is strai ghtfo rward to determine. Some of the details of the shape of these curves 
will be confirmed by calculations in section VIII. For most of the rest of the paper, we will not have much more to say 
about the phase diagram in the ($, T) plane, and refer the reader to ref. |l7j for discussions of its features^ Note, 



1 Discussed in ref. |17|, for example is the issue of the line of extremal black holes for T=0 and $>1. The calculation of 
W[T, <3>] yields a non-zero result on this line, which is the contribution from the extremal black holes. We expect that this 
does not represent the equilibrium situation, because they will decay due to "super-radiance" effects on the approach to zero 
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however, that the boundary in this figure marks the line where the Gibbs free energy of the black holes equals that 
of AdS. That is the boundary does not denote a curve where one of the local stability criteria begins to be violated. 
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FIG. 12. Sketches of the thermally stable phases in the canonical ensemble, and in the grand canonical ensemble, respectively. 



Depending upon the situation, there may or may not be the possibility of electrical fluctuations. This depends 
very much upon the setting within which we are considering these black holes. In a theory without charged particles, 
the black hole charge would be fixed and electrical stability need not be considered. In general, however, if there are 
fundamental charged quanta in the theory, then there is the possibility of the black holes emitting or absorbing such 
quanta, introducing the possibility of electrical fluctuations. Such a possibility must be considered in (for example) 
the case when the EMadS system is considered to be a Kaluza-Klein truncation of some higher dimensional theory, 
as discussed in our previous work flT]]. Then, the electrically charged black hole can in principle emit or absorb 
electrically charged Kaluza-Klein particles in order to allow its charge to fluctuate. 

In the particular case of four dimensions, however, there is also the possibility that we can exchange, by electric- 
magnetic duality, the electric charge (and vector potential) that we have been considering here for a magnetic charge 
(and vector potential). In this case, we have instead that the only way for the magnetically charged black holes to 
change their charge is to emit or absorb Kaluza-Klein monopoles, which are not fundamental quanta, as they are 
very massive, the further we are below the Kaluza-Klein scale. 

In general, when there are allowed electrical fluctuations (by whatever mechanism is appropriate to the situation 
in hand) we must also take into account on the phase diagram, the electrical stability of the solutions as determined 
in the previous section. Including those regions, we obtain the following phase diagram: 



temperature, as the charge in them is not fixed in this ensemble. This is an artifact of the failure of the Euclidean quantum 
gravity techniques that we have used to take into account such processes. 
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FIG. 13. The phase diagram in the canonical ensemble, showing the disallowed (shaded) regions where the solutions are 
unstable to electrical fluctuations. Note that the critical point and part of the coexistence line lies within the unstable region. 



The question arises as to what the equilibrium system is which resides in the shaded regions. The electrically 
unstable black holes cannot reside there, and so we must search for other possibilities. One formal possibility is that 
extremal black holes reside there, because formally they can exist at any temperature for any charge. However, we do 
not find this possibility very attractive. We expect that the permission that the Euclidean computation appears to give 
them to exist at any temperature is an artifact, and that they should naturally be associated with zero temperature, 
in which case they can only occupy the line T=0 on our phase diagram in the canonical ensemble, which they do. 
In any event, one can infer from the calculations of ref. [l7| that extremal black holes actually have a higher free 
energy than the unstable nonextremal black holes. Another possibility is that the preferred state is simply anti-de 
Sitter space (which can also exist at arbitrary temperature) filled with a charged gas. This is certainly the case at 
Q=0 |ll| , p^| . However, when the gas carries a non-negligible charge (and hence mass), its backreaction on the AdS 
geometry can not be neglected in determining the free energy. Another interesting possibility is that of black hole 
surrounded by a gas of particles. Again, if the gas component carries a sizable fraction of the charge and mass, its 
backreaction on the geometry would modify the equation of state and may then re-establish thermodynamic stability. 
Pursuing either of these possibilities lies beyond the scope of the present paper, and so we will leave settling of this 
interesting issue to a future date. Hence we must simply regard the shaded region as a sort of terra incognita with 
regard to black hole physics. As a final note, we remind the reader that this is only the region in which we are certain 
that the black holes do not minimize the free energy due to its thermodynamic instability. It may be that the onset 
of a phase transition to a state of lower free energy actually occurs outside of the boundary in figure 13, just as it 
does for the grand canonical ensemble in figure 12. 



VIII. COEXISTENCE OF PHASES AND THE CLAPEYRON EQUATION 



Let us study further the coexistence lines which we discovered in the phase diagrams, both in the canonical ensemble 
and in the grand canonical ensemble (see figure |l2| ). We can use straightforward thermodynamics to determine the 
shape of the lines. 

Let us start with the grand canonical ensemble, with Gibbs potential W[$>, T), with dW=— SdT— Qd<&. An equation 
can be derived for the line separating two phases A and B in the (<fr,T) phase diagram as follows. Along such a 
coexistence line the phases (for given ($,T)) have the same W, and so the slope of the curve Q(T) is related to the 
change in entropy and Q by: 

_ S A ~ S B 

dT Qa - Qb 1 ' 

In the case at hand, one of the phases is AdS which has zero entropy and zero charge. So we find that (for all n): 
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Here, q(T, $) is obtained from the equation of state for the corresponding branch. Equation (J3J) is the precise 
analogue of the Clapeyron equation. From it, we see that the slope of the curve is negative. For the case of n=3 we 
can give explicit expressions. In the rescaled units, we have: 



d<P _ Q 1 IT IT' 2 1 



(35) 



We see from here that the curve intersects the axes orthogonally, and its convexity, sketched in the figure [L2l follows 
from the fact that d 2 <P/dT 2 <0. 

Next (assuming the issue of electrical stability can be ignored), we consider the canonical ensemble, defined by the 
Helmholtz thermodynamic potential F[Q,T], with dF=SdT—$dQ. Along any line of coexistence of two phases, we 
have: 



dQ _ Sa — Sb 
dT ~ $ A - $B 



(36) 



The phase diagram is sketched in figure p2[ 

The Clapeyron equation can be used to find the slope of the curve at Q=0 and Q=Q C rit for the line separating the 
two black holes phases (we show the expressions for all n): 
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(37) 
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(38) 



(where we have used here that ?"™n) ~ q near q = 0, as is easy to obtain). On the scale at which we have sketched the 
coexistence curve in the previous section, it is essentially a straight line, and we have drawn it as such in figure O. 



IX. FLUCTUATIONS FOR CHARGED ADS BLACK HOLES 



In section |V|, we discussed and computed the thermodynamic quantities (specific heats and permitivity) which 
signal the stability (or not) of a black hole against fluctuations. While these quantities pertain to the response of the 
system to macroscopic thermodynamic processes which may be performed, in Euclidean Quantum Gravity, where we 
ordinarily do not have a description of the microscopic degrees of freedom, we usually cannot relate them directly to 
microscopic fluctuations, as we can in ordinary thermal physics. 

However, we can go further in this paper. Many of the adS models which we have here can be embedded into a full 
theory of quantum gravity — string and/or M-theory — and where the holographic duality tells us precisely that the 
microscopic description is organized neatly in terms of a dual (gauge) field theory 

So we may go and boldly study the fluctuations of the thermodynamic quantities in our theory, and we should see 
earmarks of the underlying (gauge) theory in our quantities, connecting the microscopic to the macroscopic. 

Here, one uses the entropy to define a probability distribution on the space of independent thermodynamic quan- 
tites pj] ]: p(Xi)ocexp(S(Xi)). With the assumption that fluctuations are small, we can work with a quadratic 
expansion of the entropy in deviations from the equilibrium values. The stability analysis of section VI establishes 
that the Hessian of S is negative semi-definite, and so we have a normalizable Gaussian distribution within this 
approximation. One then finds that the fluctuations are given by 
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where SXi denotes the deviation of Xi from its equilibrium value, and notation of the left-hand side denotes a matrix 
inverse. 

Implicit above is the assumption that we have a closed system can be divided into a number of subsystems. In 
the AdS context, the natural decomposition is the black hole and the thermodynamic reservoirs^. In this situation 
where the subsystem of interest is really the entire object under study, the most reasonable approach is to consider 
fluctuations in only the extensive variables that are free to vary in the thermodynamic ensemble [pOfl . Hence we denote 
the general extensive variables that arc free to vary as Xi, and the F^s are their conjugate intensive variables defined 



by dS=F k dX k . Eqn. (|39D then becomes j30|,|31 



<»■«"> (40) 

Now, these fluctuations are given practical meaning when they are compared to for example, their equilibrium 
values. For example, the relative root mean square of the fluctuations: 



« (41) 

which tells us about the sharpness of the distribution in Xi. 

Note, that by the formula (^) the above ratio goes roughly as the extensive parameters to the —1/2, and therefore 
the distribution is increasingly sharp as the size of the system increases. 

Now in the present problem of charged black holes, 

dS = (l/T)dE- {&/T)dQ . (42) 

Hence for the canonical (fixed Q) ensemble (our analogue of a fixed volume system), the only free extensive variable 
is the energy, and the above formulae yield 

dE\ m9 fdE\ 



For the grand canonical (fixed $) ensemble (analogue of a fixed pressure system), the energy and the charge are free 
to vary, and one has 

{SE 2 ) = - (—) = T 2 I— 

= TOs, + T$ f ) . (44) 
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(SE 5Q) = T 
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,2 i'9Q\ . _ fdQ i _ .. r . 



dT)® ' T$ \9^)t = + T$£T ' (46) 
We have recovered the fact that the thermodynamic fluctuations are controlled by the same generalized compressibil- 



ities — specific heats permittivity, etc., — that determine the intrinsic stability in section VI . This follows since both 
analyses can be phrased in terms of the Hessian of the entropy. 

Above we have presented a general thermodynamic discussion. Let us now focus on the case n=3 and present our 
results in terms of the dimensionless variables introduced in section Q. Note that translating the above thermodynamic 
formulae to the dimensionless variables, there are extra factors, giving e.g., 



lx We are neglecting the contributions of any gas component around the black hole in all of our calculations in this paper. 
Further we should be able to consider smaller subdivisions with the dual field theory in mind. 
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For the fixed charge ensemble: 



(6E 2 ) _ 3G 1 
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E 2 2vr/ 2 r\ (r\ - 1 + $ 2 ) 



(47) 



For the fixed potential ensemble: 

(SE 2 ) 3G 1 (r 2 +l-$ 2 ) 



(48) 
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(51) 



Notice that all of these results are proportional to G/l 2 ^N~ 3 / 2 , so for large N the fluctuations are suppressed. 
For n=3, the dual field theory (supplying our microscopic description) is the field theory of ref. p2| , associated to N 
coincident M2-branes. The number of degrees of freedom in this theory grows as iV 3 / 2 (as seen for example in the 
black hole entropy at high temperature). So the squared fluctuations are controlled by the inverse of the number of 
degrees of freedom of the field theory, which is precisely what we expect from standard kinetic theory connecting the 
microscopic to the macroscopic! Note here that we see these unconfined degrees of freedom appearing in our formulae 
at arbitrary temperature in this ensemble. This is because black holes dominate the thermodynamics for all values 
of the temperature: the presence of charge effects a dcconfincment of the theory at all temperatures, even in finite 
volume. (This is to be contrasted to the case of Q=0, where AdS dominates the physics for some T<Thp, representing 
the "confined" phase.) 

To gain more insight into these results, let us rewrite eqn. (|7p for the energy fluctuations as: 
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(52) 



From this form, one can pick out some of the interesting behaviour. The fluctuations go to zero at zero temperature 
as T 3 . For large T (and hence large E), the fluctuations also go to zero now as 1/T (since for large r + , 2T~r + ). An 
interesting factor is (dT/dr + )^ 1 which can change sign for Q<Q cr i t =l/(2v / 3). 

So for Q>Q cr it, the fluctuations rise from zero at the extremal black (T—0) go through a maximum and then die 
down for large temperatures. As Q approaches Q C rit, the maximum grows larger and larger, and actually becomes a 
divergence at <2=<2crit! (We have plotted these squared fluctuations in figure [l4|, where they are denoted f(T).) This 
is actually the same divergence as that in Cq at the critical point, which was commented on in ref Jl7[ |. Hence one 
finds from there that near the critical point, 
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(53) 
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FIG. 14. The squared fluctuations f(T) in the energy, relative to the equilibrium energy, for varying Q>Q cr it- The values ol 
Q plotted here are (bottom up) Q=0.49, 0.44, 0.39, 0.34, Q CYA . The dotted curve shows that the fluctuations diverge at Q—Q C iit, 
at the critical temperature T=T cr i t . 



This divergence of the energy fluctuations signals the breakdown of the Gaussian approximation considered in these 
calculations. It is also the classic behaviour of a second order phase transition point, where correlation lengths, etc., 
diverge as an order parameter vanishes. Here, the order parameter can be taken to be a homogenous function of 
r+(3)— the difference between horizon radii of the branches 3 and 1. 

For Q<Q C rit, the fluctuations rise from zero at the extremal black hole and diverge at the first zero of dT/dr+. 
Between the two zeros of dT/dr + , (5E 2 )/E 2 is negative. This is simply an indication that we are in the thermally 
unstable regime, otherwise known as "branch 2". For T larger than the second zero of dT/dr + , the fluctuations are 
monotonically decreasing (from infinity at the zero, to zero as T — > oo). As we know from the minimum free energy 
condition, we are protected from the unstable regime by making a phase transition from branch 1 and 3. So in a 
(SE 2 )/E 2 versus T plot, the fluctuations rise from zero to the phase transition point and the (discontinuously) jump 
to the decreasing curve. 



X. OF HIGHER DIMENSIONS AND OTHER THERMODYNAMIC FUNCTIONS AND ENSEMBLES 



In this section, we collect together some results for various thermodynamic quantities computed for arbitrary n, 
with all of the factors explicitly included. The thermodynamic functions are written in terms of their canonical state 
variables. We do not use the physical charge Q instead of q, for simplicity of presentation. In any expression, Q may 
be restored by recalling that 

Q = ^(n-l)cq and c = J^-^- . (54) 
87rG V n — 1 

Similarly, we also introduce the parameter s as 

Notice that it does make sense to write physical quantities in terms of q and s, since they are related to the charge 
and entropy densities. This follows from the fact that we may replace Z n-1 a> ra _i by the field theory's volume V n -i. 
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The equation of state, following from eqn. ([!(]) is: 



T = 



(n - 2)/ 2 (l - c 2 $ 2 )(c$)~ + nq~< 



Equation of state for extremal black holes is: 

2 n — 2 



q n - 



*V*e - l)(c$ e )' 



for T arbitrary. 



(56) 



(57) 



The Gibbs thermodynamic potential for the grand canonical ensemble is [17]] : 



W[T,<f>] = 



WttGI 2 



(58) 



where Q=Q(T, $) is obtained from the equation of state. Notice that W[T, $] vanishes for anti-de Sitter spacetime 
(which has Q=0), and so AdS may be thought of as the reference background for the calculation of the action, and 
indeed W^T, <&] was computed in this way in ref. Jl7| , using the background subtraction method, which we see (in the 
present work) gives the same result as the intrinsic definition by "counterterm subtraction" methods. 

The Helmholtz free energy which is the Legendre transform of VK[T, <£>] may be computed with an explicit action 
calculation, using the counterterm subtraction method to give an intrinsic ( "backgroundless" ) definition. The result 



F[T,Q] 



16ttG/ 2 



! 2 ^ 



Cf>J 



{2n-3)l 2 cq$ 



(59) 



Again, $=$(T, Q) is obtained from the equation of state. Notice that AdS with with non-zero charge is not a 
solution of the equations of motion and so cannot be considered as the "ground state" or reference background for 
this result. Indeed, this result cannot be obtained by an action calculation which uses a matching to a background, 
precisely for this reason. The counterterm subtraction technique is therefore necessary here to supply the honest 
action computation for this thermodynamic potential. It is satisfying to note that W[T, $] and F[T, Q] are Legendre 
transforms of each other, W=F+Q§, as they should be. 

We can arrive at a variety of other ensembles, with their corresponding associated potentials, by formal Legendre 
transforms. For example, we can consider the enthalpy H $], a function of entropy and potential (this notation is 
not to be confused with the Hamiltonian!). Starting from W[T, $] we can construct H—W+TS, finding 



H[S, $] = 



IGirGP 



-+l 2 



; (l-c 2 $ 2 ) 



(60) 



Notice that this function can not be obtained by performing a proper background subtraction in Euclidean gravity, 
since for any given solution we cannot find another regular solution with the same values of the entropy and the 
potential. However, the enthalpy vanishes for AdS, which could therefore be regarded as the ground state or reference 
background here. 

Another thermodynamic function in terms of its canonical variables is the (internal) energy E=W+TS+Q§, 



E[S,Q] = 



(n - l)lL> ra -l 

16ttG2 2 



(61) 



which vanishes as well for AdS. This function would define the thermodynamic potential for the microcanonical 
ensemble, and as for the enthalpy above, a calculation from Euclidean gravity should proceed by fixing the entropy 
of the state-the black hole area, if we neglect the entropy of the charged gas in AdSQ. 



2 See [B3| for work on defining the microcanonical ensemble in gravity. 
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XI. THE UNIVERSAL NEIGHBOURHOOD OF THE CRITICAL POINT 



In the section IX, we saw that fluctuations diverge as we approach the critical point in the canonical ensemble. This 



point represents a second order phase transition, as can be seen from the fact that the free energy's first derivative 
ceases to have a discontinuity there (see figures ^ and || for visual confirmation), while the divergences of the last 
section signals a discontinuity in the second derivative. 

While much of the detailed discussion of the paper has been for n=3, we emphasize here again that the results 
extend to all n>2. This is most clearly seen from the important features of the equation of state. Let us examine 
some of these more closely. 

Consider equation (|i"o|). Originating as the condition for Euclidean regularity, and hence thermodynamic equilib- 
rium, the qualitative features of /3(r + ) for varying q are shown in figure [l5| These features are the same for all n: 
There is a critical charge, q C mt, below which there are three solutions for r + for a range of values of /3, corresponding 
to the small (branch 1), branch 2, and large (branch 3) black holes, in the language of ref. [0, and in this paper. 




FIG. 15. A family of isocharge curves for the (/3,r+) form of the equation of state. Note that the middle curve is for the 
critical value of the charge, Q C rit, below which multiple branches of r+ solutions appear. The neighbourhood of the critical 
point is a universal cubic, true for all dimensions. 



That this shape persists for arbitrary n can be seen as follows. First, note that for large r + , (3(r + ) goes as ~l/r+. 
Secondly, note that the denominator of the right hand side of eqn. (|l^) after choosing scalings similar to those done 
for n=3 at the beginning of section III , is[^| 



„2n-2 



„2n- 







(62) 



which has a single positive root, r e , where p diverges. This corresponds to the T— situation, and r e is the radius of 
the corresponding extremal black hole. Given the above, any turning points for finite r + must come in pairs, and the 
condition dp/dr + —0 shows that there are only two real, positive such solutions, which we call and ^+(3), labeling 



where branch 1 ends and, respectively, where branch 3 begins. Branch 2 lies between these roots, 
determining those roots also have an elegant form (for the same rescaling as before): 



2n-2 _ 2n-4 

' _l_ ' _l_ 



+ (2n-3)g 2 = 



The equations 



(63) 



The two roots coalesce at the critical point (i.e., d 2 (3/dr+=0 also) where q=q cl -it- The value of the radius 
of this critical black hole is r + ( CI . it ) and it is at (inverse) temperature /3 C rit- For example, in the case n—3, 



3 That is, we absorb a factor of Gl 1 ^Jn/(n— 2) into Q and I 1 ^Jn/(n—2) into r+, etc. 
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the quantities {g c riti P+(crit)j Asrit} take the values {1/VT2, 1/V2, 3/(2\/2)}, while for n=4, they have the values 
{2/^/135, •v/(2/3),5/(4v6)}- The basic point here is that while the critical values themselves vary, the important 
structures do not depend upon n in any essential way. 

The neighbourhood of the critical point is extremely interesting. Because of the fact that for all n, there at most 
two turning points below q C rit, it is clear that this neighbourhood can be better written as a cubic, in terms of local 
coordinates near the point. To this end, write p=r + — r_|_( cr j t ) , P = P—Pcrit, and q=q—q cvl t, and rewrite the equation of 
state in these coordinates. The neighbourhood of the critical point is found by taking these coordinates (p, {3, q) to 
be small. 

For the example of n=3, after some algebra, we obtain 

0=fV2-^P 3 + 4p 2 + ^P- 4= + £|. (64) 



P J r /T /3 r 2V3 3/3 

Note that the quadratic and linear terms vanish with at the approach to the critical temperature (3— >0, and the term 
which contains p 3 does not vanish in this way, and so we neglect higher powers of p in favour of this one in order 
to study the near-critical behaviour. Here, and in what follows, we will also neglect terms which are not linear in q 
and p. This cubic form ( |6~4| ) may always be obtained in this limit for all n, because of the observations made in the 
preceding few paragraphs. From this, certain universal behaviour can be easily deduced, such as the critical exponent 
characterizing how fast our "order parameter" , p, vanishes. (Recall that p represents the difference in equilibrium 
radius between the black holes of branch 1 and that of branch 3; it measures the distance from the analogue of the 
"fluid" phase in liquid-gas language, where the two forms are indistinguishable.) Setting (3=0, we see that the critical 
exponent is |, since 

l) V~(Q-Qcrit)* ■ (65) 

Performing this computation for other n will change the numerical prefactor, but not the exponent, which in this 
sense deserves to be called universal. 

That a cubic equation controls the phase structure can be traced back a step further. First, notice that the three 
dimensional plot of the curve in (/3,r + ,Q) space is the cusp catastrophe, as drawn in figure [II], with some sample 
state space trajectories. We can remove the quadratic term in our cubic polynomial by shifting p by an appropriate 
amount. Multiplying overall by a normalizing factor, our cubic may be written as: 

= p 3 + Ap + B , (66) 
with A ~ 4^(3 and B ~ -0 - J~q . 

Equation ( 53 ) is actually telling us about the location of the turning points of a quartic function 

V(p) = \p 4 + ^p 2 + Bp, (67) 

where we have discarded an arbitrary additive constant. Treated as a potential, (for reasons which will be clear below), 
it is the generic form of V(p) as A and B vary that controls much of the critical behaviour in the neighbourhood of 
the critical point. (As A and B are functions of Q and T, this critical behaviour in (A, B) space translates directly 
into the earlier discussed critical behaviour in (Q,T) space.) 

The function V(p) deserves to be treated an effective potential which organizes the description of much of the local 
physics. In particular, away from the critical point, where A and B are both non-zero, the potential gcncrically 
has two minima and one maximum, the location of which are given by the solutions of our universal cubic. These 



locations may be smoothly visualized in the form of the cusp, sketched in figure 11. The location of the minima in p 
are the values ^+(1) and ^+(3), of the equilibrium black hole radii of branch 1 and branch 3, while the location of the 
maximum is ^+(2)1 the branch 2 black hole radius. The thermal stability of the branches correlates with the whether 
the turning point is a maximum or a minimum of V(p), further justifying its treatment as a potential. 

The boundary of the region where there are three solutions marks the situation where one of the minima of the 
potential V{p) merges with the maximum and disappears. This boundary is simply given by the values of A and B 
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where the cubic's discriminant, A—27B 2 +4A 3 , vanishes. (This can only happen for A<0, therefore telling us that we 
have the distinct branches below /3 C rit-) The interior of this region may be translated into (Q,T) space, where it gives 
the shaded region in the third diagram of figure ^ where branch 2 resides. 

Along the line in the (q, $) plane (or the (Q, T) plane) where B vanishes, given by q=-J (32/27)/?, the two minima 
of the potential V(p) are degenerate. This is the point at which there is phase transition, as the system moves from 
one minimum of the potential to the other. 

At the critical point, (^4=0, B=0) 7 the maximum and the two minima merge into a single minimum of the potential. 
Notice that the well formed by the potential is very flat there, and so the range of allowed fluctuations within it is 
larger at this point than at any other point in the plane, as they are less confined. We have seen this physics before 
as the divergence of the fluctuations of the microscopic degrees of freedom at the critical point. The potential V(p) 
is an effective potential for the uncharged microscopic degrees of freedom of the theory in the neighbourhood of the 
critical point. (See figure 16 for a summary of these critical points of the potential.) 




FIG. 16. The behaviour of the Landau-Ginzburg As potential at various points in the (A, B) plane. This plane maps to 
the (Q,T) plane of the charged black holes system. The line (A<0,-B=0), maps to the critical coexistence line found in that 
system. 



Also, in this language, the meaning of the swallowtail shape for the thermodynamic potential F[Q,T] is now clear: 
It is simply the actual value of the potential V(p; /?, q) at its maxima and minima: the critical line is the place where 
these two values at the minima are equal, the place where V has degenerate minima. 

This function V(p) is the A3 Landau-Ginzburg potential. The effective Landau-Ginzburg theory which we can 
write here is an effective theory of the uncharged microscopic degrees of freedom underlying the thermodynamics. 
Kinetic terms to complete the Landau-Ginzburg model would have their origins in the holographically dual field 
theory. One can in principle derive additional potential terms governing the charged degrees of freedom as well, in 
order to model the stability structure uncovered in section [v|, but we will not do that here. 

In the language of catastrophe theory |Q, the term p 4 is the basic "germ" of the cusp catastrophe, and A and B 
are the "unfolding parameters" which deform the potential, giving a line of first order phase transition points along 
the line (B=0, A<0) where its mimima are degenerate. The (A-D-E) classification of such potentials is isomorphic to 
that of certain geometrical singularities . We cannot help but wonder if this story marks the beginning of a richer 
tale involving a more profound underlying geometrical structure into which this physics is all embedded. As all of the 
physics of this paper is intimately connected to the physics of branes, perhaps the possibility of such a connection 
should be pursued. 
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